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Abstract 

We show that the tunnel coupling between semiclassical states localized in 
different minima of a smooth random potential increases when magnetic field 
is applied. This increase originates from the difference in gauge factors which 
electron wave functions belonging to different electron "lakes" acquire in the 
presence of the field. We illustrate the increase of coupling by a model calcula- 
tion of tunneling through a saddle point separating two adjacent lakes. In the 
common case, when the barrier between two lakes is much narrower than their 
size, the characteristic magnetic field is determined by the area of the lakes, 
and thus may be quite small. The effect of the field on coupling constants 
leads to a negative magnetoresistance in low-temperature conduction. 

PACS Numbers: 73.50.Jt,71.55.Jv 
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Deep in the insulating regime, the low-temperature conductivity of a two-dimensional 
electron gas is dominated by phonon-assisted electron hops between localized states. The 
conventional picture [[]J of hopping implies that localized states are formed by individual 
impurities. The lower is the temperature, the larger is the separation of the localized states 
between which a typical hop takes place. In this picture, the magnetic field affects hopping 
transport in two ways: by shrinking the wave functions of the localized states, and by 
changing the contributions of different tunneling paths to the amplitude of a hop. The 
latter means that in the course of tunneling between initial and final states an electron can 
"visit" a sequence of virtual states localized on neighboring impurities, so that the amplitude 
of a hop represents a sum of partial amplitudes. In a magnetic field, each path acquires a 
phase factor. These factors suppress the destructive interference of different paths, thus 
causing the decrease JJ,[3| of the resistance in a weak magnetic field. It was shown in that 
this decrease comes from specific hops for which tunneling amplitudes are close to zero in the 
absence of a field. Although the portion of these hops is small, the rise of conductivity with 
magnetic field is related to the increase of the tunneling probabilities for these particular 
hops, since they are most sensitive to the field. 

A different realization of the insulating regime emerges if the random potential is smooth. 
In the latter case, the electron gas breaks up into separate lakes, each lake accommodating 
many electrons. Within a certain temperature range the electron transport is provided by 
tunneling of electrons through saddle points of barriers separating the adjacent lakes. This 
picture is different from the standard picture of the variable-range-hopping |]J, so that the 
straightforward application of the theory of magnetoresistance 0|| is impossible. In the 
present paper, we study magnetoresistance of electron gas in this regime. 

It is important to note that the states active in transport (i.e., with energies close to 
the Fermi level) correspond to high-number levels of size quantization in the minima of 
the random potential. As a result, the amplitude of a hop is determined not only by the 
transmission coefficient at the saddle point, but also depends significantly on the overlap 
between the wave functions on the both sides of the barrier. This overlap is normally small 
because of the oscillatory behavior of the wave functions. If the applied magnetic field is 
weak enough, its prime role is to affect this overlap. We argue that the overlap increases 
with magnetic field, i. e., the coupling between the states, localized in adjacent minima 
of the potential, becomes stronger. In other words, unlike the conventional model, 
magnetoresistance of each elementary link appears to be negative thus leading to the net 
negative magnetoresistance of the system. 

Consider two minima separated by a barrier. In the vicinity of a saddle point the barrier 
potential has the form: 

U(x , y) = Uo -^ + ^, (1) 

see Fig. 1. We define the coupling constant between the states p and k in the neighboring 
minima as a matrix element t p k in the tunneling Hamiltonian: 

H = £ E l p a\a p + £ E r k b\b k + ^(t pfe aj6 fe + t; k b{a p ), (2) 

P k pk 

where E l , E r k are the energy levels in the two minima, and a) , b\ are the corresponding 
creation operators. The expression for t p k in terms of the wave functions ip l p {x, y) and ipl(x, y) 
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can be obtained by generalizing to the two-dimensional case of the procedure employed in 
the derivation of the tight-binding approximation ||: 




(3) 



Here m is the electron mass, and Ge{vi, IJ2) is the Green function describing the propagation 
of an electron under the barrier between the points (— d/2, y{) and (d/2, 1/2), see Fig. 1. In the 
derivation of @ we have neglected the curvature of the lines of turning points, U(x, y) = E, 
on the both sides of the barrier. The form of the barrier ([]]) enables one to separate the 
variables in the region \x\ < d/2. As a result the Green function can be expressed in terms 
of eigenfunctions (f n {y) of the harmonic oscillator with frequency u y , 

r / ^ = v ( Pn(yi)^fn(y2) U) 

E m,V2) V K(d/2)v' n (-d/2) - v' n (d/2)u' n (-d/2) ' 1 ] 

where u n [x) and v n (x) of the one-dimensional Schroedinger equation with the barrier po- 
tential Uq — muo 2 x 2 /2 and energy E — (n + l/2)%u y . If the transmission coefficient of 
the barrier is small, one can use the quasiclassical expressions for the functions in 
u n , v n oc exp{± J X dx'K n (x')}/ ^/k^ with 

1/2 

(5) 



rrvJix 1 . . 1. 

U 1 E + huj y (n + -) 



The precise form of the potential confining the electrons in each lake is a complicated 
functional of the distribution of positive charges outside the plane of 2D electron system. If 
the size of a lake exceeds the effective Bohr radius, as — fr e/me 2 , screening by the electrons 
of the lake renormalizes considerably the bare confining potential (here e is the dielectric 
constant). However, for a two-dimensional electron gas, the screened potential has still the 
spatial scale of variations equal to the lake size rather than being equal to a small screening 
radius ||. Therefore, we use the simplest parabolic form for the confinement potentials in 
order to carry out the calculations to the end. We assume 

ViM = rnnUx + D 1 + d/2) 2 + ^ (fi) 
mQ 2 rx (x-D 2 -d/2) 2 mQ 2 ry y 2 



V r (x, y) = sev ^ — + 
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for the potentials forming the left and right lakes respectively. The corresponding equipoten- 
tial lines are elliptic with excentricity depending on the corresponding ratios Q(i,r)x/Q{i,r)y 
Then the eigenfunctions in each lake are the products of the harmonic oscillator wave func- 
tions in x and y-directions. Denote with Yi and Y r the turning points of the wave functions 
in the ^-direction for the states p and k respectively (see Figure 1). If Y t and Y r exceed 
significantly the effective size of the tunneling region, the ^/-dependence of ip l p ( y —d/2,yi) and 
■0£(d/2, y 2 ) takes the simple form: 

(~^y^j K sin (j- Q iy Y iVi + Xi) , (7) 
d \ (m 



ip r k [-,y 2 j oc sin ( — Q ry Y r y 2 + Xr J , 
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where the phases of xi an d Xr are or 7r/2 depending on the parity of the numbers of 
the harmonic oscillator wave function. Obviously, such a simple form applies when both 
numbers are high. 

As an example, we discuss the coupling between two even (in the ^-direction) states 
Xi = Xr = 7r /2. Then Eq. reduces to the evaluation of the following double integral: 

t pk oc J dy 1 J dy 2 cos (^-Vti y Yiy^j cos (^-VL ry Y r y^j 

X ^2 ( Pn(yi)(Pn(y2) exp (- / dx'n n (x') \ . (8) 

V J—d/2 I 



Note that the integrals over y x and y 2 can be considered as the Fourier transforms of the 
functions (p n . Since for a harmonic oscillator the Fourier transform of the eigenfunction has 
the same functional form as ip n , we have 



tpk oc Vy?2n \ Yi— J f2n (V r — J exp ( - / dx' K, 2n {x') J . (9) 



The summation over n can be carried out if we assume the term (n + \)fiu y in Eq. (Q) to 
be a small correction and expand « n with respect to this term. Then we have 

dx'K n {x')&^ L+n( n +-)-». 10 

-d/2 Aia;^ 2 

Now we can make use of the addition rule for the functions (p n to obtain: 



t pk oc exp 



n(U - E) m(Y?nl + Y^Vl? 



ly ' r ryJ 



Tiui x 2 Tiuiy t anh ( nuj y /uj x 



( mYiY r Qi y Q ry \ 

cosh i • u / \ ■ 11 

smn(7ro; 2y /a; x .j / 



The first term in the exponent describes the transmission amplitude for the barrier at the 
central crossection y = 0. The other terms describe the reduction of t^p due to the oscillations 
of the wave functions ip p and ip^. The expansion of K n (and, hence, EqflTTD) is justified if 
the main contribution to lnt^ comes from the tunneling exponent. 

The magnetic field B affects the coupling matrix element tk p by modifying the wave 
functions ip p and ipk in the lakes and by changing the Green function Eq. (^) for the motion 
under the barrier. However, if we choose the guage in the form 

A x = 0, A y = Bx, (12) 

then the Green function will not be affected (if the width of the tunneling region is small 
compared to the lake's size). Therefore, we have to find the variation with magnetic field of 
the wave functions tp p , ipk only. Let us start from the left lake centered at x — —d/2 — D\. It 
is convenient first to find the function tp p (x, y) in the gauge A x = 0, A y = B(x+d/2+Di) and 
then to make a transformation to the initial guage Eq .([l2D. As a result of this transformation, 
the wavefunction ip p acquires a gauge phase factor: 



^ P {x,y) = 4> p (x,y)exp 



2niB (d \ 



(13) 
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$0 being the flux quantum. A similar consideration is applicable to the wave functions in 
the right lake, and those acquire an additional factor exp[— 2-niBy[dj2 + D 2 )/&o\- 

It is convenient to write the Schroedinger equation for the function ip p using the variable 
X\ = x + Di + | instead of x. Then we have 



2m 



+ 



<9 2 tL 



dx 2 dy 



— ihu c Xi 



m 



d A _ 
dy 2 



(14) 



where u c = eB/cm is the cyclotron frequency. We search for the quasiclassical solution of 
Eq. (|TJD in the form exp[iS(x,y)], where S is the action. For weak enough magnetic fields, 
uj c <C Qi x ,Qiy we can present S(x,y) as So + Si, where So is the action for motion in the 
potential m(Q 2 x x 2 + Qf y y 2 )/2 defined by the equations: 



lx 



dS mVt 



ly 



dx\ h v dy h 

Assuming Si <C S we get the following equation for the correction Sx(x,y): 
h 2 fdS dS 1 dSodS^ 



-muStiyXiJYf - y 2 . 



m \dxi dx\ dy dy 
It is easy to check that the solution of this equation has the form 

-muj c Vliy 



Si(x!,y) 



mi - n& 



VliyXxy + Vt lx JD\ 



* 2 JY 2 



(15) 



(16) 



(17) 



Eq. ( p~T| ) gives the phase which a state p acquires in the magnetic field. To calculate the 
coupling coefficient with the help of Eqs. (H) andflTSD, we need the value of Si at x x = D\. 
It is convenient to present the phase factors [originating from flTB| ) and flT7D] for the wave 
functions ip p (—d/2,yi) and 7/^(^/2,7/2) as exp(27T77/iLi/$o) and exp(— 2^7/2-^2/^0)) where 



Li, 



2 + W 



(l,r)x 



o 2 



D lt2 . 



(18) 



After introducing these factors into Eq @, the calculation of the coupling coefficient is quite 
similar to the case of zero magnetic field, and yields: 



t kp (B) = t kp (0) exp 



B 1 



2B 2 



cosh — . 



(19) 



of zero-point motion, X s = 
the saddle point of the barrier and the lakes, 



It is convenient to express the parameters B and B\ in terms of L\, L 2 , and two amplitudes 

Ti/mujy and A^ r = Jh/mtt(i <r ) x , that characterize respectively 

2LiL 2 



Bn 



$0 



+ 



\f2TiXJLl + L\ [tanh(nujy/uj x ) (L 2 + L|) smh(nujy/uj x ) 



-1/2 



(20) 



$0 



2t\L\Ij2 



A, 



Y 



— — 1 + — cosh 



LO 1 



Yr 



- ^ 1 + =± cosh 



A, 



ITU, 



Lx 



-1 



• ■. nLJ y 
smh — - . 

u)~. 
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Correspondingly the magnetoresistance between two lakes behaves as 



i2(0) V **p(°) / V B o / cosh 2 fi/Bi ' 

The magnetoresistance at small B is quadratic 

5R(B) 




(22) 



and it is negative if B > B\. To estimate the ratio B /Bi, we notice that typically all the 
dimensionless factors in (|20|) are of the order of unity, and therefore B /Bi ~ L 12 /A s . The 
latter ratio can be estimated as the number of occupied levels in each of the lakes, and is 
large in the semiclassical regime. Hence, in the small field region the resistance decreases as 
SR(B) / R(0) ~ —B 2 /B\, the characteristic value B x being of the order of & /DiD 2 . The 
resistance falls off exponentially with B in the region B c ^> B ^> Bi, where the crossover 
field is B c = Bq/Bi. After reaching a minimum at B « B c , the resistance rises sharply at 
higher fields. At B 3> B c the exponential increase of the resistance persists. The origin 
of magnetoresistance in this limit can be understood in terms of the edge states formation. 
Indeed, it can be easily seen that under the condition B 3> B c the cyclotron radius becomes 
smaller than the size of the lake. Two edge states on the opposite sides of the barrier 
correspond to the orbits skipping in opposite directions. Then the mismatch of the wave 
vectors characterizing these orbits becomes only bigger as the magnetic field increases. Note 
however that Eq is not valid at B ^> B c since at B = B c the cyclotron frequency oj c 
becomes of the order of Q x 

The physical meaning of the negative magnetoresistance is that the phase factors acquired 
by the wavefunctions in magnetic field compensate effectively the difference in their wave 
numbers in the ^/-direction and, therefore, lead to the increase of the coupling. This is clearly 
seen from the Eq. (H). Indeed, the Green function Ge(vi, IJ2) connects most effectively 
points with y x = 1/2, and is a slow function of its argument. On the other hand, 

functions ip l k (d/2,y) and if)p(d/2,y) rapidly oscillate with different periods. We have shown 
explicitly on a simple model, that a weak magnetic field reduces the "distance" between the 
Fourier harmonics belonging to these two wavefunctions thus increasing the overlap between 
dipi/dx and dip r /dx at the barrier. A similar phenomenon for the "vertical" tunneling 
between parallel two-dimensional electron systems was considered recently by Lian Zheng 
and MacDonald [0]. 

The general case of lakes of an arbitrary shape does not allow one to separe the variables, 
which may lead to "quantum chaos" ISHlOf. It might be argued [JTTJ] then, that the states 



ipi and ip r are well represented at each point by a random superposition of plane waves 
travelling in all directions. This superposition, along with other harmonics, would contain 
the plane wave striking the inter-lake barrier at right angle. Due to the latter harmonic, the 
coupling matrix elements tk iP would loose the exponential sensitivity to the magnetic field. 
However the results of the recent studies (see the book by Gutzwiller ||, and Heller's chapter 
in [12] and references therein) suggest that even in the case of chaotic motion, a substantial 
portion of quantum eigenstates with high numbers are associated with the periodic classical 
trajectories in the lake. Namely, these states can be conceived as the result of the Bohr- 
Sommerfeld quantization of plane waves propagating along these trajectories. The numbers 
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of states of the both types described above are of the same order jn| . It means that strong 
negative magnetoresistance and exponential suppression of conductivity persists for random 
shapes of the electronic lakes, because of the existence of "scarred" P,|i~2|| eigenstates. 

The magnetoresistance fl2ll) has no temperature dependence. This is because while 
deriving our main result fl20|), (|2l| ) we assumed the motion of electrons within each lake 
to be completely coherent. At a finite temperature, the phase-breaking length, L^, may 
become smaller than ■ This will eliminate effectively, the parts of the phase factors 
(|18|) acquired due to the electron motion within each of the lakes, i.e. the terms in L 12 
proportional to D 12 . Therefore at < D 12 L X 2 in ( [T%1 ) should be replaced by the value of 
d/2. As a result, in the range of temperatures where ~ Di )2 , the magnetoresistance must 
exhibit a strong temperature dependence, and become weaker at higher temperatures. The 
sign of the magnetoresistance remains negative, because for a semiclassical electron motion 
d exceeds A/ ir . jS , and the condition B 3> B 1 remains valid, even if < D 12 . 

In conclusion we have demonstrated that the hopping conductivity of two-dimensional 
electron gas subjected to a smooth random potential should increase drastically with in- 
creasing magnetic field. The larger are the electron lakes (the closer is the Fermi level to 
the classical percolation threshold), the sharper is the raise of the conductivity with mag- 
netic field. The giant negative hopping magnetoresistamce in weak magnetic fields that may 
be attributed to the mechanism we considered here, was recently reported for GaAs gated 
heterostructures [0,0. However, the observed [|l4||i"5f transition to the metallic (Quantum 



Hall) state at higher fields can not be accounted for within this simple model. 

Discussions with I. L. Aleiner, B. L. Altshuler and H. U. Baranger are gratefully ac- 
knowledged. The work at the University of Minnesota was supported by NSF Grant DMR- 
9423244. 
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FIGURES 

FIG. 1. Two electronic lakes separated by a saddle point (a), and the schematic potential profile 
in the crossection y=0 (b) 
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